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ABSTRACT 


A statistical analysis of space guidance systems, based on the 
properties of multivariate Gaussian distributions and linear perturba- 
tion theory, is discussed. The results of this analysis can be applied 
to any mission regardless of the nature or complexity of the trajectory 
or the guidance system, as long as the linear perturbation theory is 
valid. The theory needed to statistically describe the injection errors, 
trajectory dispersions at the target due to the injection errors, and the 
calculation of the average midcourse maneuver of a guided spacecraft are 
reviewed. 
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AERO-ASTRODYNAMICS INTERNAL NOTE 68-2 


THE USE OF STATISTICS IN GUIDANCE ANALYSIS 


SUMMARY 


A statistical analysis of space guidance systems, based on the pro- 
perties of multivariate Gaussian distributions and linear perturbation 
theory, is discussed. The results of this analysis can be applied to 
any mission regardless of the nature or complexity of the trajectory or 
the guidance system, as long as the linear perturbation theory is valid. 
The theory needed to statistically describe the injection errors, tra- 
jectory dispersions at the target due to the injection errors, and the 
calculation of the average midcourse maneuver of a guided spacecraft are 
reviewed . 


SECTION I. INTRODUCTION 

The application of statistics to guidance analysis has received 
considerable attention during the past few years. Noton [l] has pointed 
out that the treatment of injection errors as independent random vari- 
ables may sometimes be a poor approximation, and cross-correlations 
between injection errors must be taken into account. This paper reviews 
the theory (taking into account the cross-correlations between injection 
errors) needed to statistically describe the dispersions at the target 
and to determine the magnitude of the average midcourse maneuver for a 
typical interplanetary flight. To accomplish this, a statistical know- 
ledge of the system errors, e.g., gyro drift, accelerometer errors, 
engine shut-down errors, etc., is used to statistically describe the 
injection errors. This representation of the injection errors, along 
with a state transition matrix, is used to describe the dispersion at 
the target planet. Such information is directly related to the proba- 
bility of mission success. 

The analysis used is credited to Noton [l] and is based on the pro- 
perties of multivariate Gaussian distributions and linear perturbation 
theory. The results of this analysis can be applied to any mission 
regardless of the nature or complexity of the trajectory or the guidance 
system, as long as linear perturbation theory is valid. 



SECTION II. JUSTIFICATION FOR USING LINEAR COMBINATIONS OF SOURCE ERRORS 
TO APPROXIMATE TRAJECTORY ERRORS 


Injection error analysis is often based on the method of linearized 
perturbations. The approach taken is to represent the injection errors 
as linear functions of the source errors. A statistical description of 
the source errors will then yield a statistical representation of the 
injection errors. 

It is assumed that the injection coordinates (x, y, z, x, y, z) 
are functions of the source errors e^l g i £ n) and assume nominal 
values when e^ = 0(1 g i & n). That is, 


x « x(e, , e s , . . . , e ) 


y = y(e x , e 2 > ••• > e n > 


z ~ z ( e^ , e 2 , • • v , e ^ ) 


x - x(0, 0, . . . , 0) 
y = y (0 , 0, . . . , 0) 


z = z(0, 0, . . . , 0) 


are nominal values of the injection coordinates. 


Consider a Taylor series expansion of the first of equations (1) 
about the point (0, 0, 0). Then, 


x( e i , e 2} . 


, e ) 
’ n 


- x(0, 0, . . . , 0) + 


e i • 


• • + Se n e n + R ’ 
0 


2 


t 


• &*. jil_d k ^ 


where 


5x 

do 


0 


denotes the partial derivative of x with respect to evaluated at 
(0, 0, 0), and R is the remainder term in the Taylor series expan- 

s ion . 

For sufficiently small source errors, denoted by Se^, one can write 
since R approaches zero, 


Sx = x(e r , e 2 , e n ) - x(0, 0, 0) 


“ + a iaS e ,2 + ••* + a in § e n > 


where 


a 


c3e. 


0 


Thus, for sufficiently small e^l *5 j £ n), the injection error, 8x, is 
accurately approximated by a linear combination of the e^. Similar results 
can be shown for y, z, x, y, z. 


Then 


Sx - a^Be^. + a 12 Be 2 + . . * 4- a in oe n 
Sy = a 2 i§ e i + a 22 Be 2 + ... + a^ Se n 


h 


if : i 




5z - + a sa Se 2 + . . . + a 6n 6e n 


or in matrix notation 


SX = ABe, 


( 2 ) 


3 


where 



8x 


8®1 

8X * 

By 

• 

• 

, A * [a i{ ], Be * 

Se 2 

* 

» 


8z 


• 

Be 




n 

*•>* y» 


The terms in the A matrix would be computed on a digital computer in a 
comprehensive trajectory program. For n error sources the calculation 
of the A matrix would require n *4* .1 machine runs. On the first of these 
runs, all source errors are zero and, of course, x, y, z, x, y, z are 
nominal. On the second and subsequent runs, each independent error 
source is perturbed, one at a time, and perturbed values of x, y, z, 
y, z are found. Each such run yields a column of A. For instance, 


a 


ij 


x(0, 0, 80 ., 0, ... 0) - x(0, 0, 0) 

8e . 


a t’j 


z(0, 0, ... , Be., 0, ...0) - z (0 , 0, ... , 0) 

8e . 

J 


Thus, under the assumption that linear perturbation theory is valid, tra 
jectory errors can be accurately approximated by linear combinations of 
source errors. 


SECTION III. DERIVATION OF THE PROBABILITY DENSITY 
FUNCTION INJECTION ERRORS 


The injection errors are described statistically by their joint 
probability density function, f(SX) (see Appendix A). The derivation of 
f(f»X) is based on the properties of multivariate normal distributions 
(Appendix B ) and the results of the linearization procedure of Section II. 

The hardware errors 0^(1 ;£ i < n) are assumed to be independent random 
variables, each normally distributed with mean = 0, and variance, a?. 

The zero mean assumption is valid, since if the random hardware errors 


4 


I t /’- 1 ' i 


havr non-wro means, Lhoy conic! be corrected before launch. i« practice, 
l lie ftMiL’lKLU'H of the individual errors can be satisfactorily approxi- 
mated by norma! dls ts Lhutiona , Furthermore, since the number of system 
errors is largo and the errors arc often independent of one another, the 
central limit theorem asserts that the sum of these errors approaches a 
normal distribution in the limit regardless of the distribution of the 
individual system errors. 

Using the fact that linear combinations of independent normally dis- 
tributed random variables are normally distributed (see Appendix C), and 
the results of Section II, it follows that the joint distribution of the 
injection errors is six-dimensional normal* That is, 


-1/2(SX-B) T M" 1 (8X-B) 


where 


*2, *> ' 


E(5x) 

E(Sy) 

* 

m 

E(az) 


x-:-.. • . ■ 


: 


and H, the covariance matrix of injection errors, is given by 


M = E [ (SX) (5X) £ ] - AA A^ , (Appendix I), 


where 


0 1 0 , . . 0 


A » E[(Se)(SeX] * |0 a 2 


0 0 


A may be determined from pre-flight laboratory tests. Once A is found, 
the probability density function of injection errors can be readily 
obtained. 
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SECTION IV. DISPERSIONS AT THE TARGET 


In the absence of post injection guidance, the dispersion at the 
target is a function of the errors at injection. The statistical des- 
cription of injection errors, therefore, can be employed to obtain a 
convenient pre-flight description of trajectory dispersion at the target. 

The dispersion at the target may be measured i 'n terms of miss coor- 
dinates. The miss coordinates, M x and M I? , are measured in a plane per- 
pendicular to the approach velocity vector. 

From perturbations on the standard trajectory a matrix, U, may be 
formed such that 


M X 

Mjp 


= UbX = UASc, 


where A and de are obt ,~ned from Section II. 

Since linear combinations of independent normally distributed random 
variables are normally distributed, it follows that the joint distribution 
of M x and Mg is bivariate normal. That is, 


- 1/2 


£<Mi, M ;2 ) = 


« 1 


Mm 


^-Mm *> c U ,- 



27 C (| c |) 1 / 2 


( 3 ) 


where Mm 1 ~ Mm.., ~ 0, and C, the covariance matrix of the miss components 
is given by 

C * E[(USX) (usx) T ] * umu t , 

where M is the covariance matrix of injection errors. 
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In order to determine a convenient description of trajectory disper- 
sion, let £(Mi» M;>) a CONSTANT, Then the exponent in (3) is also a con- 
stant, That is, rewriting the exponent of (3) gives 



where lc a is a constant, p is the coefficient of correlation between Mi 
and M->, and , o M ^ are the standard deviations of M x and M a , respectively. 

Equation (4) is the equation of an ellipse (Appendix IV) having semi- 
major and semiminor axes kAi and kAs» where 


(5) 


( 6 ) 



Each of these ellipses has its major axis inclined at an angle 0 to the Mi 
axis , where 


0 


1 Can -1 

2 4 - 4 


for 4 X * a M s 


0 - it/ 4 for o' 


_ ,-2 

'Mi ' °M a ' 


(7) 


The probability of the miss being within such an ellipse is 


P - 1 - e 


-1/2 k a 


(when k = 1, 2, 3, P. « 0,40, 0.86, 0.99) (Appendix. IV) . 





SECTION V. DETERMINATION OF THE MAGNITUDE OF 
THE AVERAGE MIDCOURSE MANEUVER 


For a given mission, the analysis of the preceding section may indi- 
cate u high probability of excessively large trajectory errors at the 
target. It may then become necessary to allot propellant for making mid- 
course corrections. Let P be the probability that enough fuel is carried 
on the spacecraft to perform the required corrective maneuver. Since the 
magnitude of the midcourse maneuver depends on the covariance matrix of 
injection errors, preflight knowledge of the probable magnitude of the 
midcourse maneuver can be obtained. This makes it possible to allot a 
sufficient amount of propellant tc insure a probability, P, of accomplish- 
ing the midcourse maneuver without carrying an excessive amount of fuel. 

For a midcourse point on a given trajectory, it can be shown that, 
to a first order of approximation, the three components of the maneuver 
are given by 



V = 



= KSX = KAoe, 


where K is a (3 x 6) matrix consisting of elements computed on a standard 
trajectory. Since E(V X ) “ E(V y ) = E(V Z ) = 0, we have 


E [V V T ] = KE [BX SX T ] K T 
= KMK T . 


KMK T is the covariance matrix of the three-dimensional normal distribution 
of the components V x , Vy, and V 2 . Although V x , V y , and V z satisfy a joint 
normal distribution, the distribution of (V x 4* Vy + Vz) 1 ' 2 ' - V is not 
normal . 


Now, there is an orthogonal transformation which diagonalizes the 
covariance matrix 1<MK T so that the elements of the diagonal are the eigen- 
values of KMKT [5]. In this new coordinate system let V *i> V yi> and V z 


be the random variables denoting the components of the mid course maneuver. 
The eigenvalues of KMK^- are a 


x i* 


rf. , and a I; , the variances of V x , V v , 
yv z i x l ^l 5 


and V z respect 


ely. The following procedure (obtained from ref. 8) may 
be used to approximate the AV capability required to insure a probability, 
P, of accomplishing the mid course maneuver. 


Lot k x denote the largest eigenvalue of KMK T . If \fIT[ is at least an 
order of magnitude larger than both and a one-dimensional normal 

distribution of 



is assumed. We have 



e 


2 k- 


dx, 


whore n can be determined from Tables of Normal Distribution. Then, 

V = n \fk7- 

If n/1^i and \fk£ are of the same order of magnitude and are at least an 
order of magnitude larger than ^k 3 , a two-dimensional normal distribution 
of Vx Is assumed. Thus, in polar coordinates, 


and 


n \fk7 2rt 


P = 


1 rf 

2 k-, 


2 jtk 


rdQdr 



n = n/ ^ n(l/l-P) 2 and V = n \/kx. 


*•- r. 4*4; 


* \ * 


O 



m ■ 

■v 

t V 

t* 

rv 

i, 

7 


ir -Jk 2 , and Vk 3 are all of the same order of magnitude, a three- 

drmens lonal normal distribution of Vl is assumed. Then, in spherical 
coordinates, » AU b P n ericai 


'V/* 

if ; ; 

n \/k 7 

jc 2 jt 

V Y ' 
**. '■» * 

f 

• C 
cf ' 

p. r 

f f 1 

j V. 

0 o 

1 0 J ( 2 rit 1 ) 3 /s 

/ 

#.’• 
V ! 

n s/kT 

- I 

1 : 

•V " 

= / 

~ — e 2 dr 

£ 

-n Vkl 

4J, Tk x 


1 r 


s 


e 2 v .2 


^ sin 0 d© dfl dr 


1 dr - n e‘ n /2 


and 


V = n Vk 7 . 


#. 

®Fl 


llius , a preflight statistical analysis can be used to obtain a good esti 
mate of the magnitude of the midcourse maneuver. 
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APPENDIX A 

REVIEW OF ELEMENTARY STATISTICAL 'THEORY 

A random variable x is one whose, value can be predicted only on a 
probabilistic basis. A function f(x) > sucb that 

b 

P[x is in (a ,b)] = J f(x) dx, 

a 

is called the probability density function of the random variable x. 
The distribution function, F(x),is then 


x 


F(x) 


f(x) dx. 


A physical interpretation is given f(x) and F(x) by considering a unit 
mass distributed along a straight line such that the fraction of mass 
concentrated to the left of X = x is F(x). Then 

the density of the unit mass at point x. 

The mean or expected value of x, E(x), or p, is given by 


00 

n 


p = E(x) = 


xf(x) dx. 


-00 


E(x), the coordinate of the center of gravity of the mass distribution, 
provides a measure of the location of the distribution. 

The variance of x, V(x) is given by 

00 

V (x) = f [x - E(x)f f (x) dx; 


-po 


and again V(x) is the moment of inertia (or the second moment about the 
mean) of the associated mass distribution. V(x) is a measure of 
variation of the distribution. 

In many cases, the result of a random process is not expressed by 
one observed quantity, but by a certain number of simultaneously observed 
quantities. For example, in a study of guidance errors, the random 
experiment involves the simultaneous observation of six random variables - 
the errors in three components of position and three components of velocity. 

11 




Jpte* . -HC" • -* 1* 

* * « ' >4> 


" x , 

*- * > r# t 


The one-dimensional concepts are now generalized, and we are concerned 
with the probability that the random variables x 1 , x 3 , . . . ,x will 
simultaneously belong to the intervals (a^b,), (a R ,b 3 (a n , b n ), 
respectively. That is, 

P (a, x < b, , . . , , a < x < b ) 

— 1 ’ > n “ n — n 


> • * • dx^ » « * • >dx n f 


where f(x., , x 3 ,...,x^) is the joint probability density function of x, , . . . >x 
In general, the random variables will not be independent of one another. 

That is, Pfa. < x. < b,. ] will depend upon the values assumed by all the other 
random variables in the joint distribution. 

The mean or expected values, E^), of the x i are defined by 

t» CO 


ITK = E (x^ ) ~~ ^ ^ X ^ > • • ' ^ dxj , , , dx^ , 


-co -do 


If the second order moments are computed about the mean values, then 

the covariance a is obtained' that is 

x i x j 

CO 00 

a =E[(x. - m.)(x. - m.)] = / ... / (x. - m. )(x. - m.) 

XiX^ LV i i j j J J J l i J J 

-00 -00 

f(*i ,...,x n ) dx, ...dx n . 


If i = j , this second order moment is called the Variance c . The 

, x i 

standard deviation is simply the square root of the variance, and is 
used as a measure of the dispersion about the mean. 

Correlation coefficients of the variables x^ and x. are defined as 
fo 1 lows : 


- o - / 

x .x. x -j x . / o a 

1 J J x i x j 


The correlation coefficients are bounded between -1 and 1 and are equal 
to 1 for i = j . 
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' V * * .. • , 

"* V MkA-’ 


* V* ' >’ 


If x^ and Xj are independent, then P Xi = 0. If p ■ ±1, there is a 

complete linear dependence between the variables x^ and xj. Correlation 
coefficients provide a measure of the degree of linear dependence between 
the respective variables. 


[mxn J _ , 

1J 


The expected value of the matrix A, IE [A], is defined by 


$ 

% 


E[A] - tE(a..)]. 

X J 


X = . 


Then, 




dev.. ■ , ~ 
■Is i 


E(x x ) 
E(x) = : 

E(x n ) 


The covariance matrix, C, of the random vector x is defined by 


C = E [ (x - E (x) ) (x - E(x)) J. 


Then 


C [nxnl = E[x t - E(x.)) (x. - E( Xj ))] Inxn] 


• • -a . , 
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\ 


*:• , J,<- J f" 


■ * ,> 


or 


c = 


(T" 

X 


Px 


i x r, C *1 ° x c 


Px l x n ° x l ° x n 


symmetric 


a 


2 

X, 


Py V O’v 

x 2 x n x 2 x n x 


n 




It is also necessary to consider linear combinations of the random 
vector x and to determine the new covariance matrix. For simplicity, let 
K(x) - 0, then C = E[xx T ], Suppose 

y [mxl ] = H [mxn] ^ [nxl ] ^ 


* 

'% 


o ° ’ 

tM 

f 


The covariance matrix of y is 


M 


M = E [yy T ) 






= E{Hxx T H T } 




Let xx T =■ D; then 




M 


[mxml - E (hdh t ) 


y* ■ ■ 

• * 


= E 


n n 


>1 k=l 


h. . d h., 
rj jk ik 


where 


H-thjj], » - Idy], H 1 = [h ji ]. 
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PRECEDING 

APPENDIX B 
NORMAL DISTRIBUTION 


filmed. 


A fundamental role in probability theory is played by the function 
p(x), defined as follows; For any real number X, 


0(x) - 




1/X‘ka 
2 ( ~ ) 


o’ 


( 8 ) 


A random variable x is said to be normally distributed if its probability 
density 0(x) is given by (1) for any real number x. That is the probabil 
ity that a 5 x Xk b, 


P[a S x £ b] « 


b 

J 0(x) dx. 


a 


The parameters p. and a are the mean and standard deviation of x, 
respectively. 

The probability that x will fall in the interval p ± kcr is given 
by 


Let 


Then 


p+kcr 


P = 


_ I/£UM 
2 ' a ' 


p-kcr 


s/2 


dx. 


7tCf 


t 


X - u 
kcr 
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When k = 1, 2, 3, it is found from tables of the normal density function 
that P ® ,68, .95, and .99, respectively. 

We now consider the normal distribution of two random variables. 

The two-dimensional random variable (x, y) is said t<> have a bivariate 
normal distribution if its joint density function, f(x, y) , is given by 


1 


f(x,y) = 


1 


2 (l-p s ) cr 


x a x a y °y 


2 j«r x u y 1 - p 2 


(9) 


Tha t is , 


P[(x,y)eA] 


f f(x.y) 

A 


dy dx, 


where A cR^> |iy, a x , o y are the mean and variance of x and y, 

respectively, and p is the coefficient of correlation between x and y. 
Notice that f(x,y) may be written in the form 


- \ '|[x-p K y- My I'M" 1 


‘ X 'H> 

y-p. 


y J 


f(x,y) = 


n/ _ (2 tc) s ( 1 M | ) A / 2 


where 


M 


x 


p Vy 


p a x°y 


on 


M is the covariance or moment matrix of x and y. Again the probability 
that |j^ - g x g (jjj + and Py - o y . £ y S py + cr y is given by 


LL +CT IX + cr 

hx x ^y y 


P = 


f(x,y) dy dx. 


ll — cr u, -a 
h x X K y y 


18 


It can be shown that 


P[U* ” Ko x 2 x c n x + Ka x , n y - Kcr y si y sS My 


+ Kcf y ] * 1 


- e 


i. 

2 * 


when K “ 1, 2, 3, P * ,40, .86, and .99, respectively (Reference 4). 
The preceding discussion may be generalized to n-dimensions . Let 

X = (*i> x n , x n ) 


be a random vector. We shall say that x is normally distributed in n 
dimensions if its probability density function is 


- ■|(x-B) T M; 1 (x-B) 


( 2*) n/2 ^ 


n 


( 10 ) 


where 


B = E 


x- 


X' 


X 


n 


and Mfl is the covariance matrix of x 1( x 2 , x n< 

Equation (10) is completely determined once M n and B are known. Also, 


Pfx.eA] - J f(x) dx x ... dx n , 
A 


where A cr . 
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APPENDIX C 

LINEAR COMBINATIONS OF NORMAL RANDOM VARIABLES 


This appendix shows that a linear combination of independent 
normally distributed random variables is normally distributed. The proof 
is based on the properties of characteristic functions. 


Let *1, Xs, ...» x n be Independent normally distributed variables, 
the parameters of xj being jaj and cr. , The characteristic function of x 
is given by 


j 


itx. 

E(e J 


) 


IX. it- 

» tjr 


2 V . 

J 


a . t 


Let 


X = Xl + x 0 + ... + x . 
m ^ n 

The characteristic function of x is the product of the characteristic 
functions of all the x,. Thus, 


This is the characteristic function of a normal distribution with 
paramo tors 


“ “ L “j and 0 = 2, °y 


Thus, the sum of independent normally distributed variables is itself 
normally distributed. 

Since any linear function of a normal variable is itself normal, it 
follows, by the preceding paragraph, that a linear function 


a i x T a 2 x a + . . . 4- a n x n + b 


of independent normal variables is itself normal with parameters 


|i = a^x + a 2(i 2 + ... + a^ + b 


fr2 — -L. 4- a 2 /-y2 

a - a x o*!, i- . . . * a n CT n . 
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APPENDIX D 
DISPERSION ELLIPSES 


In this section it is shown that the probability that the miss coor- 
dinates (Mx, M 2 ) are contained in the ellipse 


k' 


1 

2 

*M 

1 

2p M. M. s 

M? I 

, . g. 

1 - p 2 

(j2 

L Mi 

a K x a M. 

"t _.2 

°TU 

m 2 j 


( 11 ) 


is given by 1 - e 


-1/2 k 2 


Equation (11) represents an ellipse in the (M l5 M 2 ) plane with center 
at the origin, or aim point, and axes that are not parallel to the coor- 
dinates axes as long as p 0. 

The coordinate system may be rotated so that the axes of the ellipse 
are parallel to the new coordinate axes. The linear transformation required 
is defined by 


where 


y 1 = M x cos 0 + M 2 sin 0 
y 2 = -Mi sin 0 +M 2 cos 0 , 


(Reference 4) 


tan 20 = 


2p a a. 

M„ M 
1 a 

2 2 

°M " °M 
1 2 


for a, £ ol, 
M 1 M 2 


and 


x £ 

0=7 for cl. = cl. . 
4 Mi M 2 


Carrying out this rotation gives 


p 2 

y 

i c 2 _ —i + — 2 
+ A 2 ’ 


where An and Ao are as defined in section III. 
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Now the quantity 



is distributed as an X s variable with two degrees of freedom. Thus, the 
probability, P, that (Mx, M 2 ) lies in the ellipse defined by (1) is given 
by 


k 2 


2r(|.i) 


where [j, = 1 . Thus, 


J' e" 1/2 X2 d (X 2 ) , 


0 


k £ 


P = 


1 e ”l/2 X 2 dx2 


0 


1 - e 


-1/2 k 2 


For lc = 1, 2, 3, P = .40, .86, and .99, respectively. 
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APPENDIX 15 
EXAMPLE PROBLEM 


In order to illustrate the use of the material presented in the body 
of this report, the following example problem is presented. The hardware 
error sources, and the matrices A, U, and K represent realistic but not 
necessarily exact values. 

The lcr values of six significant error sources (three each of con- 
stant gyro drift rate, and accelerometer constant bias error), and the 
matrices A, U, and K are given. From these the covariance matrix of 
injection errors, trajectory dispersion at the target in the absence of 
post injection guidance, and the probable magnitude of a midcourse 
maneuver for a typical interplanetary flight are determined. 

The platform hardware error sources are 


A = 



Constant Gyro 

Drift Rate 


GX 

ff 

0.025 

deg/hr 


GY 

0.025 

deg/hr 


GZ 

0.333 

deg/hr 


Accelerometer Bias 


BX 0.1 

0' 

x 10“ 6 

km/ sec 


BY 0.1 

CD 

1 

o 

t—i 

X 

km/sec 


BZ 0.1 

x 10“ s 

km/ sec 

the notation of Section II, 


.000625 

0 

0 

0 

0 

. 000625 

0 

0 

0 

0 

.001089 

0 

0 

0 

0 

i-* 

o 

i 

H 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

lCT 14 

0 


0 

0 

0 

0 

0 

10" 14 


25 


ft 


t 


The hardware errors are linearly related to the injection errors by the 
matrix A, where 


Thus , 


M = 


U = 


- 4.609 

-66.994 

3.662 

48066 x 10° 

-.27300 

-19.123 


- 5.785 

- 1.122 

-30.257 

96666 x 10 4 

-.45600 

- .24966 


- 4.104 

-17.717 

0.365 

35666 x 10° 

-.11933 

- 9.3346 


2.306 

-16.240 

28.133 

11833 x 10 7 

.12666 

8.46000 

• 

33.186 

- 0.133 

195.173 

66666 x 10 4 

.57000 

- .40666 


-26.826 

81.386 

0.520 -. 

42333 x 10° 

-.30333 

-20.623 


be found 

by the method described 

in Section 

II. 



The covariance matrix 

of Injection errors is 

given by M 

- aaa t . 


6.493 



Symmetric 

- 


-0.008 

1.021 






2.538 

0.039 

1.079 





-0.838 

-0.945 

-0.600 

1.759 



• 

-0.792 

6.309 

-0.124 

-5.964 

42.176 



0.618 

0.075 

1.091 

-2.598 

-0.591 

8.845 


The matrix U, relating injection 

errors to trajectory errors at the 

t, is taken to be 






-90.744 

-50.985 

- 1.314 

-1746.410 

- 99.936 

-.733 

-17.058 

28.050 

-92.698 

- 37.376 

-575.802 

-1873.930 


The covariance matrix of the miss components M x and Mo is given by 
.T 


C = UMU 


3,370,249 

-11,188,598 


-11,188,598 

38,535,805 


2b 



Using (3), (4), (5) of Section XV gives 



Ax - 6464.81 Ion 
Ao = 335.00 km 
0 = 16.23 deg. 


Thus, the probability of the miss being within an ellipse, centered at 
the aim point, having semimajor axis A 1# inclined at an angle 0 to the 
Mx axis, and semiminor axis A 2 is 0.40. 

The matrix K, relating injection errors to the components of the 
mid course maneuver, is 


-.216 

.0216 

-.0096 

-12736.55 

1196.78 

-537.22 

.0202 

.0560 

-.0030 

1150.61 

3317.27 

-179.25 

-.00916 

-.0029 

.1784 

-537.94 

-178.28 

10587.67 


The covariance matrix of the components of the mid course maneuver is 


T 

KMK 


.495 x 10“ os 
.383 x 10" 5 
.285 x 10“ 


Symmetric 

.422 x 10“ 5 

-.836 x 10‘ 6 .102 x 10" 4 


T 

The eigenvalues of KMK are 


kx = .11593 x 10- 4 
k £ = .90764 x 10“ 7 
k 3 = .77415 x 10“ 5 . 



• *f : 

* 



Since, \ I'k'i, are of the same order of magnitude and are at least an 

order of magnitude larger than a two-dimensional normal distribution 

of Vi is assumed. Then, to approximate the AV capability required to 
insure a probability of .99 of accomplishing the midcourse maneuver, let 
P - .99 in 


and 


n = n/ fin (1/1~P) 2 
AV - n •Jk^, 


The required AV capability is 10.33 m/sec. 
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